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It has long been conjectured that certain supersymmetric Yang-Mills (SYM) theories provide us 
with nonperturbative formulations of the string/M-theory. Although the supersymmetry (SUSY) 
on lattice is notoriously difficult in general, for a class of theories important for the string/M- 
'yr' theory various lattice and non-lattice methods, which enable us to study them on computers, have 

^ been proposed by now. In this talk, firstly I explain how SYM and string/M-theory are related. 

' Then I explain why the lattice SUSY is difficult in general, and how the difficulties are solved in 

^ theories related to string/M-theory. Then I review the status of the simulations. It is explained 

that some stringy effects are correctly incorporated in SYM. Furthermore, concrete values can be 
obtained from the SYM side, even when a direct calculation on the string theory side is impossible 
by the state-of-the-art techniques. We also comment on other recent developments, including the 
membrane mini-revolution in 2008 and simulation of the matrix model formulation of the string 
theory. 
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1. Introduction 

The second superstring revolution in mid 1990's made it clear that the lattice Monte-Carlo 
simulation can play crucial roles in order to understand important problems associated with nonper- 
turbative nature of the superstring/M theory, e.g. the microscopic origin of the black hole thermo- 
dynamics, the Big Bang singularity, the inflation, and the compactification of the extra dimensions. 
The key is a deep connection between certain supersymmetric gauge theories and superstring/M 
theory: it had been conjectured that the gauge theories provide the nonperturbative formulations 
of superstring/M-theory [1, 2, 3, 4, 5]. According to these conjectures, one can study the nonper- 
turbative aspects of string/M-theory by solving the gauge theories. In particular, the gauge/gravity 
duality conjecture [4, 5] relates the thermodynamics of the super Yang-Mills theory and the black 
hole. 

The first excitement had declined rather quickly, partly because the numerical techniques to 
study the supersymmetric gauge theories were rather premature. However now various tools are 
available, some interesting physics have been revealed, and the studies of other important problems 
are within reach. In this talk I briefly review basics of this subject. I start with giving a very 
elementary introduction of the gauge/gravity duality in Sec. 2. Then in Sec. 3 1 review the numerical 
approach. In Sec. 4 I review the relationship between the supersymmetric Chem-Simons-matter 
theory and M-theory on the multiple M2-brane background. In Sec. 5 I show an incomplete list of 
other attempts. 

Because of the page limit, my explanation will be rather brief, and I cannot cover all the 
important topics. I recommend the readers to read other reviews as well, e.g. [6, 7, 8]. 

2. The gauge/gravity duality 

2.1 Open string, closed string and supergravity 

In the superstring theory, point particles aie promoted to strings. Each oscillatory mode of 
the strings can be regarded as a particle, whose mass is heavier when the oscillation is stronger. 
There are two kinds of strings, a closed string and an open string (Fig. 1). The massless sector of 
the closed string describes the graviton, a scalar field called dilaton, and various tensor fields. The 
open string can have gauge degrees of freedom at the edges, and hence it naturally describes the 
gauge fields. Other fields e.g. the massless fermions and massive fields of the Planck scale mass 
also exist in the spectrum. The string perturbation theory describe the interaction of these fields. 
For consistency, (1 + 9) -dimensional spacetime (1 time dimension and 9 spatial dimensions) is 
required, and in order to relate it to our (1 + 3)-dimensional spacetime extra six spatial dimen- 
sions must become invisible, for example through the compactification. In this talk I consider the 
noncompact (1 +9) -dimensional spacetime unless otherwise stated. 

In the following I consider the type IIA and type IIB superstrings. There are two kinds of 
massless bosonic fields in these theories: NS-NS and R-R. The NS-NS fields include the graviton 
G^v, the rank-two antisymmetric field B^y and the dilaton (/). The expectation value of the dilaton 
is related to the string coupling constant gs, which characterizes the strength of the quantum stringy 
effects, by gs = (Fig. 2). The R-R fields ai^e rank-^ antisymmetric tensor fields C^^l..^^, where 
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Figure 1: A closed string (left) and an open string (right). Roughly speaking, the closed string describes 
graviton, while the open string corresponds to the gauge field. 



A: = 1,3 for IIA and k = 0,2,4 for IIB. 
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At low-energy, there is an effective description of the superstring theory in terms of the mass- 
less fields, which reproduces the scattering amplitudes calculated from string theory. This effective 
theory is the supergravity, which is a supersymmetric generalization of the Einstein gravity. 



2.2 Black brane and D-brane 

Type IIA supergravity has a black hole solution with the R-R 1-form charge, which is called 
the black 0-brane. The black 0-brane preserves half of the supersymmetry (16 supercharges). In 
addition, there are (1 + /?) -dimensional supersymmetric black-hole-like objects, the black ;7-brane, 
which couples to the electric or magnetic charge of the R-R {p + l)-form. Here p is even in IIA 
and odd in IIB. The metric of the black /7-brane is given by [9] 
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(2.1) 



where we have used the Yang-Mills coupling gYM and the size of the gauge group N in the corre- 
sponding super Yang-Mills theory for later convenience. A constant a' is the inverse of the string 
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tension, {t,y\\) is the coordinate of the {p + 1) -dimensions along which the brane is extended, U 
and Q. are the radial and angular coordinate of the transverse directions, and is the place of the 
horizon, which is related to the Hawking temperature of the black brane by 

T = - ^ ; ° (2.2) 

where dp = 2^^^''7r(^^^^)/^r((7 — p)/2). When p = 3, the coefficient in front of <i^^|_p becomes 
a constant independent of U, and hence the black brane geometry becomes a direct product of the 
five-dimensional spacetime described by {t,y\\,U) (AdS black hole) and the five sphere S^. When 
p y^3, there is no simple direct product structure. The dilaton, or equivalently the string coupling 
constant, is given by 

= {inf-^gl^ i^-^) ' . (2.3) 

At 7? = 3, it is independent of U . 

The black brane is a classical solution to the supergravity. Because the supergravity is the 
low-energy effective description of the superstring theory, the black brane must exist also in the 
superstring theory, at least in the parameter region where the supergravity gives a good approxima- 
tion. Then, is it possible to describe the black brane directly in the string perturbation theory? 

Actually there is a counterpart, which is the D-brane (or 'Dirichlet' brane) found by Polchin- 
ski and collaborators [10, 11]. The D;?-brane is a membrane-like object extending to the time 
and p spatial dimensions x° = t^x^ ■ ■ ,x^. (Again, p is even in IIA and odd for IIB, in order for 
the supersymmetry to be preserved.) While the closed string can propagate anywhere in the bulk 
(1 +9) -dimensional spacetime, the endpoints of the open string must be attached to the D-brane, 
or equivalently the Dirichlet (Neuman) boundary condition along x^^^ ■ ■ ,x^ {xp ,x^ ■ ■ ,x^) is im- 
posed for the endpoints (Fig. 3). 

In order to understand the physics of the black brane, one considers the emission or absorption 
of the graviton, which come from the closed string. On the other hand, the concept of the D-brane 
manifestly relies on the open string picture. Therefore, their relationship can be understood as the 
open string-closed string duality. This is intuitively very simple; as shown in Fig. 4, a loop of an 
open string stretching between two D-branes can also be interpreted as an exchange of a closed 
string. By using this simple picture, Polchinski showed that the D-brane has the same R-R charge 
as the black brane, and that it preserves the same supersymmetry, which is the evidence that the 
D-brane and the black brane are the same object. 

2.3 D-brane and super Yang-Mills 

Let us consider the low-energy effective theory of the D-branes sitting in parallel and open 
strings stretching between them. (If they are not sitting in parallel, supersymmetry is broken in 
general.) 

First let us consider the single D-brane (Fig. 3). Because the open strings are trapped on the 
D-brane, the low-energy theory is a (/? + 1) -dimensional field theory. At sufficiently low energy, 
the oscillatory modes decouple and only the massless modes remain. They are the gauge field 
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D-brane 



Figure 3: An open string must be attached to D-branes. A closed string can propagate anywhere in the bulk. 




Figure 4: A cartoon picture of the open string-closed string duality. A loop of an open string stretching 
between two D-branes can also be interpreted as an exchange of a closed string. 



^/i (M = Oi 1 ) ■ ■ ■ ) p)> scalar fields Xi{i = 1 , • • • , 9 — p) and their superpartner V^ce (oi = 1 , 2, • • • , 16). 
Note tiiat scalar fields describe tlie transverse fluctuation of the D-brane, which appears because 
the strings pull the D-brane. 

Next let us consider the system of two D-branes (Fig. 5). Now there are two types of the 
open strings: the one having both ends on the same brane, and the one connecting two different 
branes. The former is massless, while the latter has the mass, (the distance between D-branes) x (the 
string tension). Such a situation is realized in the U (2) Yang-Mills theory with adjoint matters, by 
identifying the (/,7)-component of the matrices with the strings connecting i-th and j-th D-branes. 
In fact, from the kinetic term of the U (2) Yang-Mills 

^ J dP+'x Tr {D^Xif = y dP+'x Tr {d^X^ - i[A^,Xi\)\ (2.4) 

around the diagonal background X,- = diag{xf\xf''') the mass term (xp' —xf'^Y ' I (^^j)i2p is ob- 
tained. Mass terms of the scalars and fermions appear from the interaction terms in the same 
manner. After a little bit more detailed argument, it turns out that the maximally supersymmetric 
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U (2) Yang-Mills in {p + 1) -dimension^, 

5 = I dP+'x Tr I ^-Fl, + 1 {D^,Xif - 1 + i^lfy^'D^^^f + W] } , (2-5) 

gives the low-energy effective action. When there are A'^ D-branes, the U {N) maximal super Yang- 
Mills appears. (Here 'maximal' means the theory with 16 supercharges. If more supercharges 
exist, fields with spin larger than one necessarily appear, and it is difficult to construct a consistent 
local quantum field theory.) 




Figure 5: The system of two parallel D-branes and open strings. Because string has the tension, the string 
connecting two different D-branes becomes massive. 



2.4 The gauge/gravity duality 

Now we have two descriptions for the D/j-branes: super Yang -Mills (open string picture) and 
the supergravity, or more precisely weakly coupled type IIA/IIB superstring theory, around the 
black p-brane (closed string picture). Can these descriptions become valid simultaneously? 

Because SYM takes into account only the open strings, one has to consider the limit where the 
closed strings decouple. Furthermore, because the SYM does not contain the oscillatory modes of 
the open strings, such modes must become heavy. Such a situation is realized when the D-branes 
are sitting very close to each other (the distance is much smaller than ^/a!) and the physics near 
the horizon (closer than ^/a!) is considered. In order for the weakly coupled closed string picture 
around the black brane background to make sense, the stringy corerction must be small. Therefore 
the curvature radius of the metric (2.1) must be large compared to ^/a' and the string coupling 
(2.3) must be small. When p = 3, these conditions ai^e equivalent to gyt^N 3> and g\f^ <C 1. In this 
parameter region, at least formally, SYM and weakly coupled IIA/IIB string theory picture become 
valid simultaneously. ^ 

'in this proceeding we will consider the Euclidean theories, because we want to put them on computer. Therefore 
when we wrote '(3+ I)-d' it should be interpreted as 4d Euclidean spacetime. The boundary condition is not necessarily 
thermal (i.e. anti-periodic b.c. for the fermion along the Euclidean time); we also consider the periodic boundary 
condition, which preserves the supersymmetry. 

^Of course whether this condition is sufficient for the weakly coupled description is a highly nontrivial issue- we 
are trying to apply a naive perturbative string theory picture to a very short distance ! 



6 



Monte Carlo approach to the string/M-theory 



Masanori Hanada 



Based on that observation, Maldacena conjectured these two theories are the same (Fig. 6). 
More precisely, the conjecture is as follows: 



Dp-branes 
in the superstring theory 



open string pictu 





osed string picture 



super Yang-Mills 
in (p+l)-d 



equivalent 



IIA/IIB superstring 
around black p-brane 



Figure 6: The gauge/gravity duality. 



• The IsLVge-N and strong 't Hooft coupling limit of SYM corresponds to the supergravity. 

• The large-A'^ and finite 't Hooft coupling corresponds to the classical string theory, i.e. a' > 
but g, = 0. 

• Finite-A/^, finite 't Hooft coupling corresponds to the full quantum string theory, a' > and 
gs > 0. 

By now most string theorists believe the first two statements are true, while the last one is still 
controversial. 

Because the supergravity is just a classical theory and hence much easier than SYM, people 
often solve the supergravity to learn about SYM. However, once the stringy correction is taken 
into account, the gravity side suddenly becomes very difficult. If we can somehow study the gauge 
theory side, we can make various predictions for the string theory. Here lattice theorists can play 
important roles. 

3. Super Yang-Mills on computer 

3.1 (No) fine tuning problem 

Even if a given lattice action becomes the target theory at tree level, it does not necessarily 
converge to the correct continuum limit at quantum level. In order to understand this point, let us 
recall why the Wilson's plaquette action works. The key is the exact symmetries; it has the exact 
gauge symmetry, discrete translation symmetry, 90 degree rotation symmetry, charge conjugation 
symmetry etc at discretized level. Therefore, radiative corrections must respect these symmetries, 
which assures the emergence of the correct continuum limit. If some of these exact symmetries 
are not kept, we might obtain a wrong continuum limit, which is not gauge invariant or translation 
invariant or whatever. The same problem is more familiar when one considers the fermion; with 
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the Wilson fermion, for example, the mass-parameter fine tuning is needed because of the lack of 
the chiral symmetry on the lattice. 

The trouble with the lattice supersymmetry comes from the simple fact that the supersym- 
metry algebra contains the infinitesimal translation (the spacetime derivative), which is broken on 
lattice by construction. Therefore it is impossible to preserve the whole supersymmetry algebra. 
(For detailed arguments, see e.g. [12, 13].) Still, however, various theories allow fine-tuning free 
formulations: 

• In 4d ^ = 1 pure SYM (i.e. a theory of gluons and gluinos, without matter fields), one can 
use the chiral symmetry to forbid SUSY breaking radiative corrections [14]. 

• In 2d extended SYM (i.e. the dimensional reductions of 4d o/K = 1 , 2 or 4), one can preserve a 
part of the supersymmetry unbroken. This exact supersymmetry, combined with other global 
symmetries, assures the correct continuum limit. It was first argued at perturbative level 
[15][16, 17, 18, 19] and then confirmed nonperturbatively by the Monte-Carlo simulation 
[20,21,22]. 

• For 3d and 4d SYM, matrix model techniques play important roles. In short, lattice-like 
structure can be embedded in big matrices, so that SYM theories are realized around special 
vacua of matrix models. In [23] and [24, 25], it has been shown that the (2+ l)-d and (3 + 1)- 
d U {k) gauge theories can be obtained around the ^-coincident fuzzy sphere background of 
SU {Nk) theory in (0+ l)d and (1 + l)-d, after taking an appropriate large-A'^ limit. (0+ l)-d 
and (1 + l)-d theories can be regularized as we have seen before. Another formulation can be 
found in [26], which showed an equivalence between (3 + l)-d SYM and a certain multiple 
fuzzy sphere background of the plane wave matrix model [27]. Among the formulations for 
(3 + l)-d SYM, the one in [24, 25] has advantage that finite gauge group can be realized. On 
the other hand, although the method in [26] is valid only for U (oo) gauge group, it preserve 
all the supersymmetry in a very elegant manner. (Note that, though four-dimensional lattice 
theories requires the fine tunings, however, the number of the fine tunings can be small and 
hence the fine tuning approach might be practical. See [28, 29].) 

• For Id theories, one does not have to worry about the radiative corrections so seriously^, 
thanks to the UV finiteness. It enables us to use a very efficient non-lattice method [31]. For 
a lattice approach, see [32]. See also [33] for Hamiltonian formulation. 

Although several interesting theories like the supersymmetric QCD do not allow any fine-tuning 
free formulation yet, most theories important for the string theory can be put on computer without 
any fine tuning. 

3.2 (No) sign problem 

Another possible obstacle for the lattice SUSY simulation is the sign problem, that is, the 
determinant (or Pfaffian) of the Dirac operator is complex. However, rather surprisingly, it can 
simply be avoided by using the phase quenched ensemble in all examples so far studied. 

^One must not be too naive; see [30]. 



8 



Monte Carlo approach to the string/M-theory 



Masanori Hanada 



• In 2d .J/' = (2,2) SYM (i.e. the dimensional reduction of 4d ^ = 1 pure YM), there is no 
sign problem in the continuum limit. Although the complex phase appears as a lattice artifact 
[34, 35], one can take the con^ect continuum limit without the complex phase by using the 
phase quenched ensemble [21, 22, 36]. 

• In Id and 2d maximal SYM at finite temperature (i.e. antiperiodic boundary condition for the 
fermion along the temporal circle), the Pfaffian is almost real positive at high temperature 
region. Even at rather low temperature where a good agreement with the gravity dual is 
observed, the Pfaffian is still almost real and positive [37, 38][39]. 

• With the periodic boundary condition, non-negligible phase fluctuation had been observed 
at A/^ > 3 and/or sufficiently large volume. Still, however, the phase quench simulation gives 
consistent results with the gravity dual. When the average sign is not very close to zero, the 
validity of the phase quenching can also be confirmed numerically, without referring to the 
gravity dual; see [40] for the detail. 

• For a class of Wess-Zumino models, the world-line formaUsm and the worm algorithm can 
be used to avoid the sign problem [41]. 

3.3 (0 + l)-d SYM and black hole thermodynamics 

Now let us see the simulation results and interpret them from string theory point of view. We 
start with the (0 + l)-d maximal SYM, which is (conjectured to be) dual to type IIA string theory 
around the black zero-brane. Because the coupling constant has the dimension of {massY , this 
theory is not conformal. However the gravity dual predicts a conformal-like behavior at large-A'^. 
For example, there are many massless modes whose correlators show a power-law decay, and this 
theory is deconfining (in the sense that the Polyakov loop has a nonzero expectation value) at any 
nonzero temperature. The Monte-Carlo study does not just reproduces such predictions but also 
provides us with further insights. 

• The theory has a flat direction = 0. Because of the supersymmetry, it is not lifted by 
the quantum correction. However, if one takes the initial condition as X, = 0, a metastable 
state appears. That is, eigenvalues of X, fluctuate around zero for a while and then run to 
infinity. The metastable state becomes more and more stable as N becomes larger. It is con- 
sistent with the stringy picture: because the scalars represent the DO-branes, the metastable 
state is the bound state of the DO-branes, which is regarded as the black zero-brane in the 
gravity picture. There is no instability in the black brane solution at classical level {N = oo), 
and the instability appears through the quantum effect (1 /N correction). 

In order to take the expectation value, one has to take N large enough, so that sufficiently 
many configurations are obtained before the metastable state collapses. 

• The thermal excitation energy E in SYM corresponds to the ADM mass of the black hole. 
At large-A'^, on the gravity side, it is evaluated as 



(A-'/3£) 




(3.1) 
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where a constant c has not yet determined analytically. Note that the 't Hooft coupling 
constant A = Sym^ has the dimension of (mass)^. 

Numerical study of this quantity had been performed in [37, 38, 42, 43]. In particular, in [42] 
the next leading term had been studied. By using ^^^^ = 1 Al{X-^^^Tf -^ + c{k-^^^T)P 
as the ansatz, we obtain p = 4.58(3) and c = —5.55(7), which is consistent with the analytic 
prediction p = 4.6. By a one -parameter fit with p = 4.6, the overall coefficient is determined 
as c = —5.58(1). This is a prediction from SYM, which should be reproduced from the string 
theory. 

• In the gauge/gravity duality, the supersymmetric Wilson loop which involves the scalar field 
can easily calculated on the gravity side [44]. In (0 + l)-d, one can consider the Wilson loop 
winding on the temporal circle (the Polyakov loop), which is given by W = ^TrPe^'^''^'^^^\ 
where X is one of the scalars Zi , • • • jXg. The prediction from the gravity side is (log ~ 
1.89r-^/5+C, where C is an unknown constant. The simulation result in the gauge theory 
is consistent with this prediction, by taking C = —4.95 for A'^ = 4 and —4.58 for A'^ = 6 [45]. 

• At \aige-N and at strong coupling, correlation functions can be calculated by using the 
Gubser-Klebanov-Polyakov-Witten relation [46] by solving the gravitational equations of 
motion. In [47, 40], various two-point functions had been calculated numerically on the 
Gauge theory side, and perfect agreement with the dual gravity prediction [48] had been 
observed akeady at A'^ = 2. 

3.4 (1 + l)-d SYM and black hole/black string transition 

Let us consider the (1 + 1) -dimensional U (N) maximal SYM on spatial with circumference 
L. We impose the periodic boundary condition both for the bosonic and fermionic fields. This 
theory describes the system of Dl-branes winding on the circle. By taking the T-duality along 
the circle, it can also describe the system of A'^ DO-branes. The compactification radius becomes 
L' = Ina' /L. The positions of the DO-branes correspond to the eigenvalues of the Wilson line 
along the spatial circle, W = Pe'/'*>'^''. When the Wilson line phases clump to a point ((W) / 0), 
they describe a bunch of DO-branes localized along the , which is the black hole (black 0-brane) 
in the gravity picture. On the other hand, when the Wilson line phases are uniformly distributed 
on the unit circle ((W) = 0), the DO-branes form the (uniform) 'black string'. (Note that this black 
string is an object in type IIA string, which is different from the black 1-brane in type IIB.) One 
can also consider other objects like nonuniform black string and multiple black hole. 

Let us fix the temperature and change the radius of the spatial circle. When the radius is 
sufficiently large, the black hole solution in the noncompact space provides us with a good approx- 
imation. As the radius becomes as small as the size of the black hole, the compactified direction 
gets completely filled by the black hole; then the black hole winds on the circle and turns into the 
black string. This is the black hole/black string transition [49]. In the gauge theory, this transition 
can be identified with the breaking of the Z^v center symmetry W — )■ e^^^/'^W , because W is an 
order parameter of the center symmetry (Fig. 7). By looking at the eigenvalue distribution of W as 
a function of temperature T and compactification radius L, one can study the stringy effect on the 
transition. 
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Figure 7: The center broken phase (left, bottom) and unbroken phase (right, bottom) correspond to the black 
hole (left, top) and uniform black string (right, top). The dots in the bottom figures represent the Wilson line 
phases, which take the values on the unit circle. 

At low temperature, the phase transition can be studied from the gravity side, by comparing 
the free energy of the black zero-brane and black string solutions to type II A supergravity [50]. 
(At small volume, strings winding on the compact circle become light, so that the IIB supergravity 
approximation is not good. Then by taking the T-dual, the winding modes are mapped to the KK 
momentum, and description with IIA supergravity becomes valid. It can be shown that the phase 
transition takes place where type IIA supergravity, rather than type IIB supergravity, gives a good 
approximation.) It turns out the black hole (center broken phase) is preferred at small L, and there 
is a first order transition to the uniform black string (center unbroken phase). 

At high temperature, because all the KK modes along the temporal direction (and hence all 
the fermionic degrees of freedom) decouple, the system reduces to the bosonic matrix quantum 
mechanics (bMQM). Note that the spatial in the two-dimensional theory coiTcsponds to the 
temporal in the bMQM; hence the deconfinement transition of the latter can be regarded as a 
remnant of the black hole/black string transition. See [50, 51] for detail studies of bMQM from 
this point of view. 

In [39], the phase structure had been studied by using the twisted SUSY formulation. Al- 
though larger scale simulation with larger A'^ and finer lattice is needed for rigid conclusion^, the 
result seems to interpolate known results in high and low temperature regions, and looks rather 
convincing. It is very important to pursue this direction further by using several lattice formula- 
tions. 

3.5 Toward the simulation of (3 + l)-d maximal SYM 

Although a simulation of the four-dimensional maximal SYM is very expensive, there is steady 
progress. 

• In [52, 53], a numerical simulation had been performed by using the large-A'^ equivalence 
between a special vacuum of a (0+ 1) -dimensional matrix model and the 4d SYM [26]. 

^Because N is not large enough to stabilize the flat direction, a small scalar mass is introduced. 
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They had evaluated certain Wilson loops [52] and correlation functions [53], and confirmed 
the prediction of the AdS/CFT duality. 

• As argued in [36], although all known lattice formulations have a fine-tuning problem, the 
lattice simulation may still be manageable. After the conference an interesting paper [29] 
appeared, in which a lattice simulation had been performed with this philosophy. As a first 
guess, the simulation has been performed without any parameter fine tunings. At all bai^e 
couplings they studied, their results look consistent with the theoretical predictions. The 
next important step should be to see the restoration of the supersymmetry, with a possible 
parameter fine tuning. 

• Four-dimensional SYM can also be regularized by combining the two-dimensional lattice 
and fuzzy sphere [24, 25]. Because the same technique is applicable for the 3d SYM [23], 
which is numerically much cheaper, it should be easier to test the validity of this method at 
nonperturbative level by a simulation of 3d SYM. 

• It has been argued that physics of the supersymmetric sector in the strong coupling region 
can be captured by a simple matrix model [54]. For numerical work along this direction, see 
e.g. [55]. 

4. ABJM and membrane mini-revolution 

The gauge/gravity duality is expected to hold also in M-theory, whose precise definition is 
yet to be known. In paiticulai^, a stack of multiple M2-branes (i.e. the usual supermembranes in 
eleven dimension [56]) should be described by a three-dimensional supersymmetric conformal field 
theory, whose dual gravitational description is the M-theory on the AdS^ x background. Such a 
theory had been found in 2008 by Aharony, Bergman, Jafferis and Maldacena (ABJM) [57]. It is a 
U {N) X U {N) Chern-Simons theory with bifundamental matter fields, which is rather complicated. 
Its coupling constant, or equivalently the Chern-Simons level, is quantized to be integer k. M- 
theory on AdS^ xS'' /T^u gives a good description when N ^ , while type IIA superstring theory 
on AdS/^ X S'^ /CPt, is good when k <^ N <^ k^. The dual M-theory geometry predicts a peculiar 
scaling of the free energy, F ~ 

It looks hopelessly difficult to study this system on lattice. However, rather unexpectedly, 
Monte Carlo approach turned out to be very effective for studying this system. 

The key is the localization method [58, 59], which maps the ABJM theory to a much simpler 
matrix model. Let us consider the paitition function. 



Suppose the action and the path-integral measure are invariant under a fermionic symmetry gen- 
erated by Q. We consider a one-parameter deformation, 5, = S + tQV, where V is an arbitrary 
fermionic function and 0>t, and consider Z, = J[d^]e^^'. Then 




(4.1) 



Zx{Q{Ve-'Q''))=0, 



(4.2) 
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and hence Zj is actually f -independent. In the limit f — )• oo, tQy gives an infinitely steep potential, 
the path-integral 'localizes' to the classical solution to this potential, and one-loop calculation 
around the classical solution provides us with the exact answer. With a clever choice of V , the 
classical solution becomes very simple, sometimes as simple as one matrix, and then the partition 
function Z is the same as that of the matrix model. In 4d maximal SYM, for example, the resultant 
matrix model is simply Gaussian [58]. The same argument can be repeated for the expectation 
values of the Q-invariant operators; for example, a supersymmetric Wilson loop in 4d maximal 
SYM can be calculated by using the Gaussian matrix model. 

The localization method had been applied to the ABJM theory in [59]. The partition function 
can be written as 



-imyliinrii-Y n,-i2cosh^^^ ' 



n,<,(2sinh^)'(2sinh^'' 



I(M?-v?) 



(4.3) 



Analytic studies of this theory confirmed F ~ \/kN^/'^ including the prefactor [60] , and impor- 
tant result on the 1/A'^ correction had been obtained [61]. On the other hand by using the standai^d 
Monte Carlo technique this integral can be studied much more easily [62]; because the simula- 
tion is so cheap, we can obtain quite accurate numerical values. As a result, by combining the 
analytic and numerical techniques the analytic expression to all order in 1 /A^-expansion, which 
corresponds the string perturbation theory, had been obtained [61, 62]. After the lattice conference, 
analytic techniques to study this integral has been developed further; exact values are obtained for 
various and k [63], and by subtracting the perturbative contribution the instanton effect has been 
discussed [64] . Now the gauge theory side is under perfect control, and the quantum coiTcction 
on the gravity side is going on, in order to test the AdS/CFT duality at fully quantum level. It is 
interesting to apply similar numerical techniques to other important theories. 



5. Other topics 

There are many other interesting developments in the numerical simulations. Here is an in- 
complete list^ : 

• 4d jV = 1 minimal super Yang-Mills has been actively studied [65]. Probably the first target 
is the spontaneous breakdown of the 'L2N symmetry (here A'^ is the size of the gauge group 
SU{N)) due to the gaugino condensation. (Interesting arguments on SYM multiplet on the 
lattice [66] appeared after the lattice conference.) 

• Spontaneous SUSY breaking is also a very important subject. Hori and Tong [67] con- 
jectured that in the two-dimensional c/K = (2,2) SYM, which is obtained from 4d ^ = 1 
minimal SYM through the dimensional reduction, the supersymmetry is spontaneously bro- 
ken. Although numerical work by Kanamori et al [68] is consistent with unbroken SUSY, it 

^Let me know if you find something missing here. 
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is not clear whether they observe the same vacuum Hori and Tong considered^, and further 
study is needed. 

• The Wess-Zumino models on lattice are also actively studied. See [69, 70] for the restora- 
tion of supersymmetry in the continuum limit, [71] for the RG flow, [72, 73] for a use of 
the Nicolai map and the relation to a certain superconformal minimal model, [74] for an 
improvement of the lattice actions of the matrix quantum mechanics, and [75, 41] for the 
spontaneous SUSY breaking. 

• In [76] the ground state structure of the IKKT matrix model [2], which is the dimensional 
reduction of the 4d ^ = 4 SYM to zero dimension, has been studied. In this model, the 
distribution of the eigenvalues may be regarded as the spacetime; there are ten matrices, and 
if the eigenvalue distribution is four-dimensional it might be regarded as the emergence of 
our four-dimensional spacetime. In a sharp contrast to previous studies, they considered the 
Lorentzian version directly, and found that one time direction and three expanding spacial 
dimensions appear. 

Emergence of the spacetime from the point of view of the noncommutative geometry is also 
actively studied. (See e.g. [77, 78] for recent progress. ) It would be interesting to apply 
numerical techniques to such scenarios. 

Acknowledgement 

I would like to thank M. Honda, I. Kanamori, J. Nishimura and S. Shiba for useful comments 
on the draft. 

References 

[1] T. Banks, W. Fischler, S. H. Shenker and L. Susskind, Phys. Rev. D 55, 5112 (1997). 

[2] N. Ishibashi, H. Kawai, Y. Kitazawa and A. Tsuchiya, Nucl. Phys. B 498, 467 (1997). 

[3] L. Motl, ai-Xiv:hep-th/9701025; R. Dijkgraaf, E. P Verlinde and H. L. Verlinde, Nucl. Phys. B 500, 43 
(1997). 

[4] J. M. Maldacena, Adv Theor. Math. Phys. 2, 231 (1998). 

[5] N. Itzhaki, J. M. Maldacena, J. Sonnenschein and S. Yankielowicz, Phys. Rev D 58, 046004 (1998). 
[6] S. Catterall, D. B. Kaplan and M. Unsal, Phys. Rept. 484, 71 (2009). 
[7] J. Nishimura, arXiv: 1205.6870 [hep-lat]. 
[8] S. Catterall, arXiv: 1005.5346 [hep-lat]. 

[9] G. W. Gibbons and K. -i. Maeda, Nucl. Phys. B 298, 741 (1988); G. T. Horowitz and A. Strominger, 
Nucl. Phys. B 360, 197 (1991). 

[10] J. Dai, R. G. Leigh and J. Polchinski, Mod. Phys. Lett. A 4, 2073 (1989). 

^This theory has several (perhaps metastable) vacua [21]; in the one Hori and Tong considered, eigenvalues of the 
scalar field are bounded around zero. 



14 



Monte Carlo approach to the string/M-theory 



Masanori Hanada 



[11] J. Polchinski, Phys. Rev. Lett. 75, 4724 (1995). 

[12] M. Kato, M. Sakamoto and H. So, JHEP 0805, 057 (2008). 

L13J G. Bergner, JHEP 1001, 024 (2010). 

[14] D. B. Kaplan, Phys. Lett. B 136, 162 (1984); G. Curci and G. Veneziano, Nucl. Phys. B 292, 555 
(1987). 

[15] D. B. Kaplan, E. Katz and M. Unsal, JHEP 0305, 037 (2003). 

[16] A. G. Cohen, D. B. Kaplan, E. Katz and M. Unsal, JHEP 0308, 024 (2003); JHEP 0312, 031 (2003). 
D. B. Kaplan and M. Unsal, JHEP 0509, 042 (2005). 

[17] R Sugino, JHEP 0401, 015 (2004); JHEP 0403, 067 (2004); JHEP 0501, 016 (2005). 

[18] S. Catterall, JHEP 0411, 006 (2004). S. Catterall, JHEP 0506, 027 (2005). 

[19] A. D'Adda, 1. Kanamori, N. Kawamoto and K. Nagata, Phys. Lett. B 633, 645 (2006). 

[20] 1. Kanamori and H. Suzuki, Nucl. Phys. B 811, 420 (2009). 

[21] M. Hanada and 1. Kanamori, Phys. Rev. D 80, 065014 (2009). 

[22] M. Hanada and I. Kanamori, JHEP 1101, 058 (2011). 

[23] J. M. Maldacena, M. M. Sheikh-Jabbari and M. Van Raamsdonk, JHEP 0301, 038 (2003). 

[24] M. Hanada, S. Matsuura and F. Sugino, Prog. Theor. Phys. 126, 597 (2011); M. Hanada, JHEP 1011, 

112(2010). 

[25] T. Takimi, JHEP 1208, 069 (2012). 

[26] T. Ishii, G. Ishiki, S. Shimasaki and A. Tsuchiya, Phys. Rev. D 78, 106001 (2008). 

[27] D. E. Berenstein, J. M. Maldacena and H. S. Nastase, JHEP 0204, 013 (2002) [hep-th/0202021]. 

[28] S. Catterall, E. Dzienkowski, J. Giedt, A. Joseph and R. Wells, JHEP 1104, 074 (2011). 

[29] S. Catterall, P. H. Damgaard, T. Degrand, R. Galvez and D. Mehta, arXiv: 1209.5285 [hep-lat]. 

[30] J. Giedt, R. Koniuk, E. Poppitz and T. Yavin, JHEP 0412, 033 (2004). 

[31] M. Hanada, J. Nishimura and S. Takeuchi, Phys. Rev. Lett. 99, 161602 (2007). 

[32] S. Catterall and T. Wiseman, JHEP 0712, 104 (2007). 

[33] J. Wosiek, Nucl. Phys. B 644, 85 (2002) [hep-th/0203116]. 

[34] J. Giedt, Nucl. Phys. B 668, 138 (2003). 

[35] H. Suzuki, JHEP 0709, 052 (2007). 

[36] S. Catterall, R. Galvez, A. Joseph and D. Mehta, JHEP 1201, 108 (2012). 

[37] K. N. Anagnostopoulos, M. Hanada, J. Nishimura and S. Takeuchi, Phys. Rev. Lett. 100, 021601 

(2008). 

[38] S. Catterall and T. Wiseman, Phys. Rev. D 78, 041502 (2008). 

[39] S. Catterall, A. Joseph and T. Wiseman, JHEP 1012, 022 (2010). 

[40] M. Hanada, J. Nishimura, Y. Sekino and T. Yoneya, JHEP 1112, 020 (2011). 



15 



Monte Carlo approach to the string/M-theory 



Masanori Hanada 



[41] D. Baumgartner and U. Wenger, arXiv: 1104.0213 [hep-lat]; D. Baumgartner, K. Steinhauer and 
U. Wenger, PoS LATTICE 2011, 253 (2011); U. Wenger, Phys. Rev. D 80, 071503 (2009). 

[42] M. Hanada, Y. Hyakutake, J. Nishimura and S. Takeuchi, Phys. Rev. Lett. 102, 191602 (2009). 

[43] S. Catterall and T. Wiseman, JHEP 1004, 077 (2010). 

[44] J. M. Maldacena, Phys. Rev. Lett. 80, 4859 (1998); S. -J. Rey and J. -T. Yee, Eur. Phys. J. C 22, 379 
(2001). 

[45] M. Hanada, A. Miwa, J. Nishimura and S. Takeuchi, Phys. Rev. Lett. 102, 181602 (2009). 

[46] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. Lett. B 428, 105 (1998); E. Witten, Adv. 
Theor. Math. Phys. 2, 253 (1998). 

[47] M. Hanada, J. Nishimura, Y. Sekino and T. Yoneya, Phys. Rev. Lett. 104, 151601 (2010). 

[48] Y. Sekino and T. Yoneya, Nucl. Phys. B 570, 174 (2000). 

[49] R. Gregory and R. Laflamme, Phys. Rev. Lett. 70, 2837 (1993). 

[50] O. Aharony, J. Marsano, S. Minwalla and T. Wiseman, Class. Quant. Grav. 21, 5169 (2004). 

[51] N. Kawahara, J. Nishimura and S. Takeuchi, JHEP 0710, 097 (2007); G. Mandal, M. Mahato and 
T. Morita, JHEP 1002, 034 (2010). 

[52] M. Honda, G. Ishiki, S. -W. Kim, J. Nishimura and A. Tsuchiya, PoS LATTICE 2010 (2010) 253; 

[53] M. Honda, G. Ishiki, J. Nishimura and A. Tsuchiya, PoS LATTICE 2011, 244 (2011). 

[54] D. Berenstein, JHEP 0601, 125 (2006). 

[55] D. Berenstein, R. Cotta and R. Leonardi, Phys. Rev. D 78, 025008 (2008). 

[56] E. Bergshoeff, E. Sezgin and P K. Townsend, Phys. Lett. B 189, 75 (1987). 

[57] O. Aharony, O. Bergman, D. L. JaEferis and J. Maldacena, JHEP 0810, 091 (2008). 

[58] V. Pestun, Commun. Math. Phys. 313, 71 (2012). 

[59] A. Kapustin, B. Willett and I. Yaakov, JHEP 1003, 089 (2010). 

[60] N. Drukker, M. Marino and P Putrov, Commun. Math. Phys. 306, 5 1 1 (201 1); JHEP 1111, 141 
(2011). 

[61] H. Fuji, S. Hirano and S. Moriyama, JHEP 1108, 001 (2011); M. Marino and P Putrov, J. Stat. Mech. 
1203, P03001 (2012). 

[62] M. Hanada, M. Honda, Y. Honma, J. Nishimura, S. Shiba and Y. Yoshida, JHEP 1205, 121 (2012). 

[63] Y. Hatsuda, S. Moriyama and K. Okuyama, JHEP 1210, 020 (2012); P Putrov and M. Yamazaki, 
Mod. Phys. Lett. A 27, 1250200 (2012). 

[64] Y. Hatsuda, S. Moriyama and K. Okuyama, arXiv: 121 1.1251 [hep-th]. 

[65] J. Giedt, R. Brower, S. Catterall, G. T. Fleming and R Vranas, Phys. Rev. D 79, 025015 (2009); 

M. G. Endres, Phys. Rev. D 79, 094503 (2009); K. Demmouche, F. Farchioni, A. Ferhng, I. Montvay, 
G. Munster, E. E. Scholz, J. Wuilloud, Eur. Phys. J. C69, 147-157 (2010); G. Bergner, I. Montvay, 
G. Munster, D. Sandbrink and U. D. Ozugurel, arXiv: 1210.7767 [hep-lat]; S. W. Kim et al. [JLQCD 
Collaboration], PoS LATTICE 2011, 069 (2011). 

[66] H. Suzuki, arXiv: 1209.2473 [hep-lat].; arXiv:1209.5155 [hep-lat]. 



16 



Monte Carlo approach to the string/M-theory 



Masanori Hanada 



[67] K. Hori and D. Tong, JHEP 0705, 079 (2007); 

[68] I. Kanamori, F. Sugino and H. Suzuki, Prog. Theor. Phys. 119, 797 (2008); I. Kanamori, Phys. Rev. D 
79, 115015 (2009). 

[69] G. Bergner, T. Kaestner, S. UMmann and A. Wipf, Annals Phys. 323, 946 (2008); 

[70] D. Kadoh and H. Suzuki, Phys. Lett. B 684, 167 (2010); Phys. Lett. B 696, 163 (2011). 

[71] T. Kastner, G. Bergner, S. UMmann, A. Wipf and C. Wozar, Phys. Rev. D 78, 095001 (2008); 
M. Heilmann, D. F. Litim, F. Synatschke-Czerwonka and A. Wipf, arXiv: 1208.5389 [hep-th]. 

[72] H. Kawai and Y. Kikukawa, Phys. Rev. D 83, 074502 (2011). 

[73] S. Kamata and H. Suzuki, Nucl. Phys. B 854, 552 (2012). 

[74] S. Schierenberg and F. Bruckmann, arXiv: 1210.5404 [hep-lat]. 

[75] C. Wozar and A. Wipf, Annals Phys. 327, 774 (2012); 

[76] S. -W. Kim, J. Nishimura and A. Tsuchiya, Phys. Rev. Lett. 108, 011601 (2012). 

[77] H. Steinacker, Class. Quant. Grav. 27, 133001 (2010). 

[78] A. Chatzistavrakidis and G. Zoupanos, SIGMA 6, 063 (2010). 



17 



